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Modelling of the magnetic accretion flow in HU Aquarii 
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ABSTRACT 

We present a magnetic stripping model for AM Her type objects. Our model is based on 
an equilibrium condition between ram pressure and magnetic pressure in a stiff dipolar 
magnetic field. We investigate the detailed geometry of the stripping process, most of 
which can be tackled analytically. By involving additional numerical calculations, the 
model allows the prediction of phase-resolved spectra and Doppler tomograms. The 
emission line features from the companion star, the horizontal stream and the accretion 
curtain are identified with the emission line components found by Gaussian fitting to 
observational data of HU Aqr in its high accretion state. Given the simplicity of the 
model, its agreement with the observation is remarkably good and enables us to derive 
a number of physical quantities such as stellar masses and radii, total mass accretion 
rate, bulk temperature, and coupling density. 

Key words: accretion - magnetic cataclysmic variables - AM Her objects - modelling 
- stars: individual: HU Aqr 



1 INTRODUCTION 

Magnetic cataclysmic variables (MCVs) are binary systems 
consisting of a late-type near main sequence secondary star 
and a magnetic white dwarf primary. Matter lost by Roche 
lobe overflow from the secondary star interacts with the 
strong magnetic field of the white dwarf before being finally 
accreted. These systems are ideal laboratories for studying 
the interaction of supersonic gas streams with strong mag- 
netic fields. 

In an AM Her object, or "polar", the magnetic field 
of the synchronously rotating white dwarf is sufficiently 
strong to prevent the formation of an accretion disk. A rela- 
tively bright representative of this class is the eclipsing polar 
HU Aqr with an orbital period of P = 0. d 08682. It has been 
under intensive study (e.g. Glenn et al. 1994, Hakala 1995, 
Schwartz et al. 1995) since its identification in 1992 (Hakala 
et al. 1993, Schwope et al. 1993). The strength and orienta- 
tion of the primary's dipolar magnetic field has been deter- 
mined with fair precision from cyclotron harmonics observed 
in low resolution spectra, from polarisation measurements, 
and from the X-ray light curve; it was found for the mag- 
netic field strength and the position of the accretion spot, 
respectively, B s w 35MG, tp s w 40°, tf s « 26°. The mass 
ratio and inclination of the system remain somewhat more 
uncertain (Q := M wd /M rd £ [2.5,5.7], i € [81°, 90°]), but 
the eclipse phase width A<jf> = 0.0763 however imposes a 
tight relationship between Q and i. Doppler maps gener- 



ated from high resolution trailed spectra by Schwope et al. 
(1997a) clearly show that the illuminated face of the sec- 
ondary star's Roche lobe and the accretion stream in the 
orbital plane are the dominant emission-line features of the 
high accretion state. The observed structure of the hori- 
zontal stream indicates that the stream initially follows the 
ballistic trajectory and we are therefore dealing with the 
situation described by Liebert & Stockman (1985), in which 
thermal pressure dominates magnetic pressure at the inner 
Lagrangian point L\. Recent observations suggest that this 
no longer holds for the reduced accretion state, when the 
density of the stream is significantly lower (Schwope et al. 
1997b). Then, as proposed by Mukai (1988), the flow may 
couple directly to the magnetosphere at the Li-point as as- 
sumed in early AM Her models (e.g. Schneider & Young 
(1980) for VV Pup). 

In this article we develop a simple but predictive model 
for phase-resolved emission line profiles and Doppler tomo- 
grams of MCVs. We test this model on observations of 
HU Aqr in its high accretion state. The accretion stream 
starts at the inner Lagrangian point L\ from the Roche lobe 
of the secondary star as shown in Figure The initial part 
of the stream, not yet affected by the magnetic field, is con- 
fined to the orbital plane and is hence called the horizontal 
stream. In our calculations it is assumed to have a Gaussian 
density profile that is however truncated when the density 
drops below a critical value below which the matter couples 
to the magnetosphere. This stripping threshold, displayed 
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Figure 1. Sketch of the HU Aqr system projected on the x-y- 
plane (top) and x-z-plane (bottom); see text for details. 



by the solid lines around the stream centre in Figure |l|, is 
defined by the ram pressure in the stream dropping below 
the magnetic pressure of a stiff dipole field. The stripped ma- 
terial, as it leaves the orbital plane, forms a two-dimensional 
accretion curtain and glides along the field lines until it is 
finally accreted on the white dwarf's surface. Even though 
magnetic stripping occurs virtually over the whole length of 
the horizontal stream, most of the matter is stripped near 
the end of the horizontal stream, as indicated by the dashed 
lines in Figure |IJ The curtain trajectory passing through the 
end point of the horizontal stream has been named vertical 
stream in the literature. The dotted line marks the trajec- 
tory of matter that is stripped right at the Li-point (see 
below) . 

X-ray radiation from the accretion region excites line 
emission from the companion star surface and the accre- 
tion stream. The line emission from the secondary star was 
found to be optically thick in our initial study (Schwope 
et al. 1997a), whereas for the emission from the horizontal 
stream and the curtain the situation is less obvious and we 
are here comparing the limiting cases of optically thin and 
thick radiation reprocessing. 

We first focus in Section ^ on the horizontal stream 
and derive the stripping rate with which matter is removed 
from it by the magnetic field. In Section ^| we investigate the 
geometry of the accretion curtain in which the matter prop- 



agates through the magnetosphere. Our model provides the 
velocity vector and densities at each point of the flow, en- 
abling us to model phase-resolved spectra and Doppler maps 
in Section which are fitted to observations of HU Aqr in 
Section |H| This allows us to derive the physical quantities 
of the system such as stellar masses and radii, total mass 
accretion rate, bulk temperature, and coupling density. Fi- 
nally in Section ^| we discuss the discrepancies of model and 
observation and propose further improvements. 



2 MAGNETIC STRIPPING 
2.1 The horizontal stream 

As suggested by the high-state observations, the stream is 
not strongly affected by the magnetic field as it passes the 
inner Lagrangian point L\. This situation invites use of the 
semi-analytical solution for an inviscid isothermal flow in a 
Roche potential given by Lubow & Shu (1975, 1976) with 
additional corrections for the magnetic influence. First the 
applicability of the Lubow-Shu-model on AM Her stars de- 
serves some discussion: (a) The stream and secondary star 
are heated by ionising radiation from the accretion region. 
The diminished line emission of the secondary's leading 
side implies considerable shielding by the accretion curtain 
(Schwope 1997a), which must therefore have a high optical 
depth for the ionising radiation. The column density there 
is even lower than in the horizontal stream, where we conse- 
quently also expect a high optical depth. This confines the 
irradiative heating to the skin layers and causes a highly in- 
homogeneous temperature distribution over the stream. We 
assume that the bulk of matter in the stream centre remains 
at a constant temperature and the overall behaviour of the 
stream is not much affected by the hot surface layer. We also 
neglect the impact of radiation pressure on the stream, (b) 
In the region where the magnetic pressure already overcomes 
the thermal pressure but not yet the ram pressure, which is 
the case for almost the whole stream apart from a very thin 
shell around the Li-point, instabilities shatter the stream 
to small blobs (Liebert & Stockman 1985). We assume that 
the mean density on a scale comparable to the stream width 
is still described by a simple Gaussian profile as originally 
proposed by Lubow & Shu. In the Reynolds averaged sense 
we treat the ensemble of gas blobs as a quasi-steady flow 
and assume a vanishing eddy-viscosity therein, (c) Further 
downstream magnetic stripping erodes the stream from the 
outside. This violates the assumption of a Gaussian density 
distribution and might also affect the temperature structure. 
Even though the stream motion is highly supersonic the in- 
ternal gas pressure is important in the transverse directions 
of the flow. Our model however assumes for simplicity that 
the magnetic pressure exerted on the surface of the stream 
is equal to the pressure that the stripped layers would have 
exerted, so that the kinematic behaviour of the core is not 
affected by the missing envelope. 

Following the suggestion of Lubow & Shu (1975) we 
formulate the problem in dimensionless units of length, time 
and mass: 

1[L] = d , 1[T] = ft' 1 , 1[M] = MO.' 1 , (1) 

where d is the orbital separation of the system, P = 2tt/Q 
its rotation period, and M the total mass accretion rate. 
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Also from Lubow & Shu (1975) we adopt the treatment 
as a perturbation analysis in a small parameter e, which 
is the isothermal sound speed yJkT/m in the stream for a 
temperature T and a mean molecular weight m in above 
units given by Equation (Q). Under the assumption of an 
isothermal flow the perturbation parameter e is constant 
over the whole accretion stream. For the orbital period of 
the HU Aqr system we have in case of a pure hydrogen gas 
e = 0.011(T/10 4 K) 1/2 (d/10 9 m)- 1 . 

Initially we consider a stream without a magnetic field. 
We use the natural coordinate system of the horizontal 
stream where s measures the distance from the Li-point 
along the stream centre and the pair (n, z) represents the 
distances from the stream centre parallel and perpendicu- 
lar to the orbital plane. Lubow & Shu (1975, 1976) propose 
that the density in the horizontal stream has a Gaussian 
distribution around the ballistic trajectory, i.e. 



P(s,n,z) = po(s) exp(- — (jn 2 + X z 2 )) 



where the characteristic horizontal and vertical widths 1/^/7 
and 1/s/Xi respectively, are functions of s only and can be 
found from an initial value problem as the solution of a 
closed set of non-linear singular differential equations, in 
which the mass fraction p appears as a parameter; the mass 
fraction is defined as the white dwarf mass over the total 
mass: p := Af wd /(M wd + M ld ). A proper numerical inte- 
gration of these equations requires asymptotic analysis near 
the singular point L\ (see Appendix |X|). With u s o as the 
s-component of the stream velocity in the centre we have in 
lowest order of e 



Po(s) 



HX 



2%e 2 u a 



(3) 



since we have normalised our equations to the total mass 
flow, i.e. J pu s o dndz = 1 and the local stream velocity u s 
can be approximated by the central velocity u s o in lowest 
order of e. 



2.2 The stripping process 

The ram pressure of the stream orthogonal to the magnetic 
field can be written as 

1 2-2 

Pram = ~ PU a0 Sill ZD 
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where zu is the stripping angle between stream motion and 
the magnetic field. Further the thermal pressure of the gas 
is given by 
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We now compare the total gas pressure p t h + Pram in the 
stream with the local magnetic pressure p mas = B 2 /(8n). We 
assume then that the stripping threshold is given by p ta + 
Pram = Pmag- Since we are considering the stripping process 
at large distances from the dipole source we approximate 
the local B by the magnetic field So at the stream centre 
with an error of order e. Using Equation (^) & (|^) we find 
that the critical surface is a centred ellipse with the principal 
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Figure 2. Stripping threshold (solid lines) in n and z direction, 
respectively for e = to 0.24 in steps of 0.04; the remaining 
parameters are as in Table |lj. Superimposed (dashed lines) arc 
stream lines for the 10% to 90% (in steps of 10%), 97.5% and 
99.9% percentiles. 



axes dn and a 2 aligned to the coordinate axes perpendicular 
to the direction of propagation: 



(6) 



with 



In 
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(7) 



For numerical calculations we assume a dipole field 



B 



3r wd (r wd -M) 



(8) 



with r wd the distance from the white dwarf centre and A4 
the magnetic moment. Figure ^| shows numerically calcu- 
lated streamlines and stripping surfaces for the geometry 
of HU Aqr. Minima in a n and a z occur immediately af- 
ter the Li point for e^0.08, i.e. for a cold stream. They 
arise because the thermal pressure decreases due to the di- 
lution caused by the acceleration of the stream before the 
ram pressure increases enough to take over the pressure bal- 
ance with Pmag that provides our stripping condition. The 
locations of the minima are close to the equilibrium point be- 
tween thermal pressure and magnetic pressure, downstream 
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from which Liebert & Stockman (1985) postulate shatter- 
ing of the stream into fine blobs. The presence of minima in 
the stripping thresholds suggests a 'magnetic nozzle' at the 
beginning of the horizontal stream, however one should be 
cautious with that interpretation since the diamagnetism 
of the secondary star, and possibly its intrinsic magnetic 
field, induce in the vicinity of L\ a large distortion of the 
ideal dipole field assumed in our model. The observed trailed 
spectra of HU Aqr show no evidence for a manifest stripping 
directly at the Li-point. Cropper (1986) however concludes 
from polarisation measurements in ST LMi that a two-spot 
accretion is the best fit to his data, so stripping directly at 
the Li-point may occur in other AM Her objects. Since we 
can fit the parameter e only very indirectly (c.f. Section |s|) 
it is very reassuring that towards the end of the horizontal 
stream, where we expect the largestproportion of the mag- 
netic stripping to occur (c.f. Figure H), the solution becomes 
virtually independent of e. 

Neglecting diffusion at the edges of the stream, the 
shape of a truncated Gaussian density profile is pre- 
served as the stream propagates. The maximum mass 
flow that can survive magnetic stripping and proceed be- 
yond position s on the horizontal stream is given by 
J(_"_)2 + (^)2 < i P u sO dndz . Consistent with neglecting 
diffusion, we consider only the mass flow m on those stream 
lines that have not yet coupled to the magnetosphere. This 
introduces a minimum over all previous stream coordinates 
s. This has a significant effect only for e^$0.08, for which 
stripping near L\ occurs. Evaluation of the integral in ellip- 
tical coordinates and substitution of Equation (^|), (^) and 
(H) results finally in 
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(9) 



The symmetry of the integration region causes the first 
order error terms to cancel. The resulting stripping rate 
A = —dmjds is then given by 



A : 



d ( e B u s0 
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(10) 



wherein Bo, u a o, 7 and \ are an regarded as functions of s. 

At the Li-point we have a mass flow m slightly below 
1 since the magnetic pressure is very small in comparison 
with the gas pressure there but it is not vanishing totally. 
This is not a major inconsistency in our model since, as dis- 
cussed above, our approximation of a dipole field is already 
very inaccurate in this region. However for sensible values 
of e the stripping occurs mainly at large distances from the 
Li-point, where our approximation holds well (see Figure^). 
The allowed mass flow m is monotonically decreasing by con- 
struction of Equation (Q) , so un-physical negative stripping 
rates are ruled out. The above calculation did not account 
for blocking of the stripped matter by the dense core of the 
stream as proposed by Mukai (1988). Surface currents in- 
duced in the stream lead to a large local distortion of the 
threading field which then voids the stream; hence our model 
assumes that matter stripped of the back side is led around 
the stream by the distorted field lines. 
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Figure 3. The stripping rate A according to Equation ( |To[ ) for 
the HU Aqr geometry (see Table Q); t takes the values 0.04, 0.12, 
and 0.36. 



3 GEOMETRY OF THE ACCRETION 
CURTAIN 

3.1 The threading process 

Liebert & Stockman (1985) point out that the threading of 
matter onto the magnetic field involves a variety of physical 
processes, which are not yet fully understood (c.f. Cropper 
1990, Warner 1995). However it seems clear that in accre- 
tion onto compact objects magnetic instabilities play the 
key role. Burnard et al. (1983) have shown that Kelvin- 
Helmholtz instabilities caused by a shear flow with inho- 
mogeneous density distribution generate small drops which 
are easily captured by the magnetic field. It is also suggested 
that, even though the gas stream is highly supersonic, the 
negative density gradients additionally produce Rayleigh- 
Taylor instabilities. By these instabilities the stream is 
broken up into large diamagnetic blobs that, in contrast, 
are only slowly penetrated by the magnetic field. We as- 
sume here that the stripping is primarily driven by Kelvin- 
Helmholtz instabilities. In our model the droplets stripped 
from the surface follow abruptly the magnetic field lines and 
their kinetic energy from the motion perpendicular to the 
field is instantaneously thermalised. It turns out that the 
resulting rise in temperature allows the flow to become sub- 
sonic again. 



3.2 The magnetic curtain 

In the previous section we have modelled the rate at which 
material from the horizontal stream is captured by the mag- 
netic field. The geometry of the threading region depends 
on the mass ratio Q and on the orientation of the magnetic 
dipole {{>d, <p<i) relative to the canonical Cartesian system of 
the co-rotating frame. In polar coordinates (r, 9, ijj) aligned 
with the dipole we use the usual parameterisation of the 
field lines in azimuth ip and equatorial radius r m ; a dipole 
field line is then described as a function of the co-latitude 8 
by 



ip = const 



(11) 
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For an arbitrary threading point (x,y, 0) on the hori- 
zontal stream we now consider the spherical triangle with 
the dipole axis, its projection onto the x-y plane, and the 
threading point as corners. The angle ip c from the projected 
dipole axis to the threading point is given implicitly by 



tan((/p<i + (p c ) = 



/' 



(12) 



wherein the quadrant is provided by the standard sign con- 
vention for the denominator. Spherical trigonometry pro- 
vides the coupling co-latitude 6 C and azimuth ip c in dipole 
coordinates 



= cos if c sin t>d 



sini/j 



(13) 



which is unique since in our case always 9 C £ [0,7r[ and ij) € 
[— -f , ■§[, it is also clear that always | sin^c] < | sin^cj. The 
above set of equations seems to us slightly easier than the 
one given by Ferrario (1989), but nevertheless describes the 
same geometry. Equation (ill]) finally yields the equatorial 
radius of the threading field line 

(14) 



sin 8 C 



We obtain the initial velocity with which the stripped 
material moves on the curtain from the purely kinematic 
argument that the component parallel to the field lines is 
conserved. For the dipole orientation in HU Aqr (see Ta- 
ble ^) it turns out that the material turns a rather sharp 
corner of zu ~ 75°. If we use a mass specific Roche potential 
normalised with an additive constant such that $(£1) = 
the dissipated energy per unity mass is given as — $ c sin 2 vj, 
where <E> C is the potential at the coupling point. Towards 
the end of the horizontal stream the energy loss for an H- 
atom is nearly 1 keV. We assume for our model that the 
resulting rise in temperature is relaxed instantaneously to 
the isothermal equilibrium temperature characterised by the 
isothermal sound speed e - in practice a different value of e 
may apply to each field line in the curtain. 

Continuity and momentum conservation for an isother- 
mal flow inside a flux tube provide a differential equation 
for the velocity u along the field lines, given the magnetic 
field B and the Roche potential $ 



du f<? \ _ e 2 dB 1 
u \ v? J u 2 B v? 



(15) 



Shortly after coupling the stream becomes highly supersonic 
again, hence pressure effects are negligible everywhere and 
we get by integration 



$) 



(16) 



which is simply energy conservation for a free falling particle. 
In principle the field lines could lead so far uphill in the 
Roche potential that Equation (|l^) predicts the material 
to stagnate somewhere with zero velocity. For the HU Aqr 
geometry this case never occurs in our simulation. 

Gliding down the field lines the material from the 
stream forms an essentially two-dimensional curtain. We 
now introduce a curtain coordinate system (c.f. Figure ^ 
which consists of the distance s of the coupling point from 
Li, the magnetic co- latitude 6, and the distance 6 perpen- 
dicular to the curtain. The coordinates are orthogonal at the 
coupling point and the curves of constant b are field lines. 




Figure 4. Sketch of the horizontal stream with a part of the 
accretion curtain, which illustrates the coordinates introduced in 
the text. 



Since the curtain flow follows field lines, we can represent 
the density as 



p = <7 c (s,0) /(A _1 6) 



(17) 



where the typical width A of the curtain depends on s and 
6. The convergence of field lines along the curtain flow can 
be expressed by the variation of the thickness of a flux tube 
perpendicular to the direction of motion. For the flux tube 
starting at stream coordinate s the distance dl to the adja- 
cent field line starting at s + ds is given by 

(dl\ 2 , 2 (dr m \ 2 2 /cW .... 

where the derivatives (drm/ds) and (dtp/ds) are a direct 
consequence of Equation (ff2|)-(|l4|) with x = x(s) and y 
y{s). The metric coefficients fo 2 .. 
tion (| 



and (B7) as functions of 9 and r„ 



and are given by Equa- 



The mass flow transported on the flux tube defined by 
s and s + ds is given by Ads (c.f. Equation (|lo|)). As a 



consequence of continuity we have then 



A (ds 
u\dl 



A 



Or 



n2 . ,2 ,9lp 2 

-5—) +hAs~) 
as v as 



(19) 



at the coupling point we know from simple geometry that 
(dl/ds)o = seen; holds. Since the flow is bound to a flux 
tube, we can write the typical width of the curtain as A -1 oc 
B(dl/ds). A quite suitable choice for / could be a Gaussian 
profile with a dispersion A, for which at the coupling point 
A 2 = a n a z holds. 



4 SIMULATION OF DOPPLER MAPS AND 
TRAILED SPECTRA 

4.1 Companion star line emission 

Our model assumes that the companion star exactly fills its 
Roche lobe, that the reprocessed radiation is proportional 
to the absorbed X-rays, and that the emission is diffuse ac- 
cording to the foreshortening given by Lambert's law, i.e. 
with a cosine dependency of the viewing angle. Isotropic re- 
emission has shown to provide a less accurate line profile. 
Thermal broadening at the bulk temperature causes a ve- 
locity dispersion e^ld/c, which is just above the numerical 
resolution of our spectral grid. The X-rays are assumed to 
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Figure 5. Model line emission from various parts of the system displayed as trailed spectra (first column) and in the Doppler tomograms 
(second column) calculated by a filtered back-projection. Rows: (a) companion Roche lobe, (b) horizontal steam, (c) accretion curtain. 
The gray-scale of each row is chosen separately. Model parameters as given in Table |lj 



originate from a point source at the white dwarf's centre. 
The involved error in the location of the X-ray source is of 
the order of the white dwarf radius i? wd , i.e. typically 1% of 
the orbital separation d. 

The upper panel of Figure [B| shows the emission arising 
of the X-ray heated face of the secondary Roche lobe. To 
match the orientation of Doppler maps usual in the litera- 
ture we invert the sign convention of the v x and v y axes with 
respect to the x and y coordinates, respectively. The numer- 
ical calculation makes use of the fact that obviously the sec- 
ondary star is totally opaque to X-rays and to optical line 
emission. The resulting self-occultation of the Roche lobe 



should provide a constraint on the orbital inclination, but 
the endpoints of the narrow emission line (NEL) eclipse are 
superimpose with emission from the stream. We therefore 
use the inclination determined for a given mass ratio from 
the eclipse width in the optical light curve. The K-velocity 
and the width of the secondary star emission feature can be 
used to estimate the relative Roche lobe size, and hence for 
the mass ratio. The maximum width of this feature is given 
as 2.45A by Schwope et al. (1997a), which unfortunately is 
not far from the spectral resolution (1.6A). 
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Figure 6. Synthetic trailed spectrum of reprocessed line emission from the secondary star in assumption of a fully transparent stream 
in X-rays and line. This is to be compared to the middle row of Figure |H| Doppler tomogram. 



4.2 Emission of the horizontal stream 

The horizontal stream emission is shown in the middle panel 
of Figure ^[ Our model assumes that the horizontal stream is 
opaque to X-rays but transparent to the optical line photons. 
For phases smaller than approximately 0.55 we are viewing 
the irradiated side of the stream facing the X-ray source; 
in contrast, for phases larger that about 0.55 we are facing 
the un-irradiated side, however the observed emission in the 
Hell 4686A line is just slightly lower by about 20% in the 
observations (c.f. Schwope 1997a) . So the horizontal stream 
cannot be opaque for both the X-ray irradiation and the op- 
tical line emission. Simulated tomograms (see Figure ^ show 
that the stream cannot be fully transparent in both spectral 
regimes either. In that case the volume emissivity would be 
proportional to the square of the local density, which pro- 
duces too much emission close to the Li-point. It has been 
observed that the Roche lobe is significantly shielded from 
the X-ray irradiation by the accretion curtain, which is even 
more dilute than the horizontal stream. So we suggest a hor- 
izontal stream opaque to the ionising radiation but nearly 
fully transparent in the Hell line so that we observe the ir- 
radiated front side even when viewing through the stream. 
The low line opacity needed seems plausible, given the large 
velocity gradients in the stream. 

These arguments make it plausible that the line radia- 
tion mainly stems from the skin layers of the stream, which 
probably is surrounded by some kind of hot X-ray heated 
corona. If we assume that the radiation reprocessing involves 
the same physical process as on the companion Roche lobe, 
then comparison of the relative fluxes allows an estimate of 
the stream width, which is proportional to e. Knowledge of 
e in turn allows an estimate of the core temperature of the 
stream. 

Gaussian fitting of the high velocity component (HVC) 
emission shows clearly a fading of the horizontal stream 
around phases and 0.5 (gaps). This suggests a foreshort- 
ened optically thick emission. However the above argument 
gives evidence for an optically thin stream transparent to op- 
tical line photons. We therefore propose a thin X-ray heated 



region that is optically thick and reprocesses most of the 
ionising radiation whereas the remainder of the stream is 
optically thin. The resulting spectra are shown in the mid- 
dle panel of Figure ^. We do not expect information from 
the gaps to be of great importance for the parameter fit- 
ting because of the rather poor signal to noise ratio at these 
phases. However, compared to the observational data they 
occur at quite the right phase and velocity in our simula- 
tions. The rather narrow emission lines (knots) at <p « 0.22 
and tp ~ 0.72 occur when the orientation of the stream to the 
line of sight is such that the intrinsic velocity plus the orbital 
speed are virtually constant over the length of the stream. In 
addition at ip w occultation by the secondary star shields 
the stream from the observer. Whereas the amplitude and 
the phases of the maxima of the horizontal stream emission 
can be adjusted to the observational data by an appropri- 
ate choice of the model parameters, the line profile of the 
high velocity wings provides a test to the correctness of our 
stripping model. 



4.3 Emission of the accretion curtain 

The broad baseline component (BBC) emitted by the ac- 
cretion curtain (see bottom panel of Figure ^) appears very 
dim in the modelled trailed spectra and in the Doppler tomo- 
grams since it is distributed over a large range of velocities. 
We assume the line emission is proportional to the inter- 
cepted X-rays. As for the horizontal stream the eclipse of 
the curtain is modelled with a spherical companion star of 
the effective Roche radius. The brightest part of the curtain 
is the stagnation region in which the material enters the cur- 
tain when it couples to the magnetosphere. Emission from 
the curtain partially covers up the gaps of the horizontal 
stream emission in the trailed spectrum of the middle row 
of Figure 0, in which the three components (Roche lobe, 
horizontal stream, and curtain) are merged. 
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Figure 7. Comparison between the observational data of HU Aqr in He II 4686 (upper row) and the model folded with the spectral 
resolution (middle row), the residuals are shown in the bottom row. Displayed are the trailed spectra (l.h.s.) and back-projected Dopplcr 
tomograms (r.h.s.), apart from the top right panel where MEM tomography was used. Model parameters from Tabic pi 



5 TESTING THE MODEL 

To adjust the 10 free parameters of the model, which are 
given in Table jj], we have performed a y^-fit to the HU Aqr 
trailed spectra from August 1993, when the object was in its 
high accretion state. We focus on the 4686A He II line which 
has the best signal to noise ratio and shows the sharpest fea- 
tures of all lines in the observed spectral range from 4180 A 
to 5070 A. Figure compares the observational data and 
our fitted model, which was folded by the observational res- 
olution. The fitting employed a simplex algorithm to adjust 
8 parameters in order to minimise the x 2 statistic for which 
we have used estimated standard deviations from a detector 



noise model. The orbital period P is very well known from 
long-term light curves; also from the optical eclipse the or- 
bital inclination i of the system is given as a function of the 
mass ratio. So we have not varied these parameters in the 
fit. 

There are certain features in the spectra which con- 
strain particular parameters. With the period and mass ra- 
tio fixed the orbital separation d scales the wavelength axis. 
The mass fraction \x governs the width and amplitude of 
the Roche lobe emission. These two parameters hence are 
predominantly determined by the NEL feature, but the cur- 
vature of the horizontal stream also effects these parameters. 
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Table 1. Model parameters obtained from x 2 -fitting; these are 
used for the figures shown in this article. The quantities where no 
error is given have not been varied in the x 2_n t- 



parameter 


symbol 


value 


error 


orbital period 


P 


0. d 086820446 




orbital inclination 


i 


84.°0 




orbital separation 


d 


5.663- 10 10 cm 


5.6-10 8 cm 


mass fraction 




0.7950 


0.0027 


magn. momentum 


eM 


0.077f 


0.0035 


dipole azimuth 




37°8 


l.°6 


dipole co-latitude 


#d 


12°4 


0.°3 


sound speed 


e 


0.0355 


0.0016 


rest wavelength 


A,) 


4685.7A 


o.iA 


rel. curtain intens. 


Ic 


0.0154 


0.0023 



The resulting Q ratio of 3.9 lies in between the two values 
obtained by Schwope (1997a) from either of these features 
separately. With d and /i fixed, the length of the horizon- 
tal stream is determined by the magnetic dipole parameters 
M, (fid and fid- However the magnetic field strength varies 
strongly with he distance from the white dwarf and only 
weakly with the dipole orientation, thus the later cannot be 
well constrained alone from a fit to the trailed spectrogram. 
Therefore we employ the knowledge of the accretion spot 
position (ifs = 40°, $ s = 26°), which is given by Schwope et 
al. (1997a), based on the X-ray bright phase and cyclotron 
harmonics. We have used the mass radius relation for white 
dwarfs given by Pringle (1975) to constrain our variation of 
ipd and "dd to those values that are in agreement with the 
observed accretion spot. The sound speed e is a measure of 
the bulk temperature, which scales the geometrical width of 
the horizontal stream. Thus its total brightness provides the 
information to which e is fitted. The model also contains the 
relative brightness I c of the curtain and the rest wavelength 
Ao as free parameters. 

The error estimates given in Table |l] have been obtained 
by bootstrapping the original observation. On each pixel a 
random signal was added with the detector noise as stan- 
dard deviation. We have repeated the x 2_n t starting with 
the same initial values for 20 different samples. The errors 
in Table [i] are ler intervals of the fitted parameters and con- 
sequently are only a measure of the uncertainty due to the 
noise in the data. They certainly do not account for system- 
atic discrepancies of the model from reality (see Section^). 

From the fitted model parameters we are able to de- 
rive a couple of physical parameters of the HU Aqr sys- 
tem. With the observed magnetic field at the accretion 
spot of 35 MG we are able to scale our intrinsic units 
from Equation (Q) to physical units. This enables us to 
determine the total mass accretion rate as M = 2.9 x 
10~ 10 M Q /yr(B/35MG) 2 (e/0.0355) 2 . This value depends 
strongly on our assumption of identical effectivity of radi- 
ation reprocessing on the companion star surface and the 
horizontal stream. Here also the choice of the stripping con- 
dition and the assumption of a quasi-steady flow are very 
essential. Nevertheless the derived M is not extremely far 
from M w 0.5 x 10~ 10 M Q /yr estimated from the total lu- 
minosity of an average CV just below the period gap (Pat- 
terson 1984). The accretion rate in the model is lowered by 
a factor of 4 when either the radiation reprocessing in the 



Table 2. Some physical quantities of HU Aqr as derived from the 
model parameters given in Table 111 



quantity 


symbol 


value 


mass ratio 


Q 


3.9 


dip angle 


</>di P 


-0.117 


primary mass 


M„ d 


O.75M 


secondary mass 


M rd 


0.19M Q 


white dwarf radius 




7.3-10 8 cm 


secondary Roche radius 


Rrd 


1.5-10 10 cm 


bulk temperature 


T 


1710 3 K 


total mass flow 


M 


2.9-lO~ lo M /yr 


coupling density 


n 


7.6 ■ 10- n g/cm 3 



horizontal stream is twice as effective as on the Roche lobe, 
or if matter coupling to the magnetosphere needs a magnetic 
field twice as large as assumed. 

In the optical light curve a pre-eclipse dip centred at 
phase (f> — —0.118 has been observed; it has a symmetric 
distribution that can well be approximated by a Gaussian 
profile with a dispersion of 0.015. Our model cannot ac- 
count for this dispersion, which is presumably a consequence 
of the matter being broken up into blobs of different sizes, 
hence the threading points vary over a range of locations. 
The good agreement of the modelled turn-off point at phase 
<f> = —0.117 shows that the horizontal stream is free falling 
in a good approximation. So the geometrical location of the 
coupling region and the fitted location in velocity space co- 
incide. 

The model also gives strong evidence that the compan- 
ion star is on the main sequence. The radius of a main se- 
quence star with a mass of O.19M0 is given by Neece (1984) 
as 0.23Rq which agrees very well with the effective radius 
of the Roche lobe provided by Eggleton (1983) which comes 
out to be 0.22Rq for our values of \i and d. We do not have 
to assume a mass radius relation for this result. 

At the coupling point we have a mass density of 7.6 x 
10 _11 g/cm 3 or, for a pure hydrogen gas, a number density 
of n — 4.2 x 10 16 cm -3 , which results in a Debye length of 
\d = 4.1 x 10 _3 cm. Hence we have an ideal plasma with a 
plasma parameter F := l/(nX d ) = 3 x 10~ 10 . 



6 DISCUSSION 

In the comparison between our model and the observational 
data of Figure ^ five details differ significantly: (1) At the 
location of the Roche lobe feature the residual tomogram 
shows a dark spot at positive v x and a bright region towards 
negative v x . This is a consequence of the missing shielding of 
X-ray flux by the curtain in our model. (2) The horizontal 
stream in the observation is clearly shifted in direction of 
smaller v y . This is not a tomographic artefact as it is also 
clear from the residuals in the trailed spectrum near the knot 
at (j) = 0.22 that the simulation and observation are shifted 
against each other. The reason for this might be a magnetic 
drag force that acts transversely on the horizontal stream 
before the matter is totally coupled to the field. (3) The 
curtain emission seems to be too faint at low velocities in 
the model spectra and also ranges out to too high velocities. 
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Also the curtain is located at lower v y in the model than in 
the tomogram of the observed data. This could be explained 
either by a reduction of radiation reprocessing when the cur- 
tain gets closer to the white dwarf, or by the importance of 
pressure effects when the curtain is funneled towards the 
accretion spot, or as well by a distortion of the dipole geom- 
etry. There also is a noticable difference in the shape of the 
eclipse, so it might be suggested that the reason lies in the 
geometry of the curtain flow. (4) The observed tomogram 
shows faint but significant emission in the region roughly 
between (v x ,v y ) = (-900,0) to (-1200,-600), which can- 
not be explained in our model. A plausible speculation on 
the origin of this emission could be larger blobs in the hor- 
izontal stream that are sufficiently massive to survive the 
stripping process but feel a strong drag force by the mag- 
netic field that forces them to leave the ballistic trajectory. 
(5) Even though folded with the observational resolution the 
widths of all modelled stream features are too small in both 
the trailed spectra and the Doppler tomograms. This could 
be caused by an underestimate of the observational resolu- 
tion, or an intrinsic line broadening. Thermal broadening 
of the bulk temperature is already included in the model, 
the gas is much too dilute for pressure broadening, and the 
field strength is much to low for a noticeable Zeeman ef- 
fect. However the layer in which the radiation reprocessing 
occurs may be heated far above the bulk temperature. This 
hot corona-like gas around the stream could be the source of 
the large observed line width. Also broadening by turbulent 
gas motion might be important. 

We have developed a predictive model for the stripping 
of the accretion stream and the line emission in AM Her 
stars in general and applied it to spectral observations of 
the eclipsing polar HU Aqr in its high accretion state. The 
main observed features, including the length of the horizon- 
tal stream and the main emission line components, are well 
reproduced by our model. A x 2_ fit of synthesised to observed 
trailed spectrograms has uncovered a number of physical pa- 
rameters of the binary system (mass flow rate, stream tem- 
perature, stream density, orientation of the magnetic field, 
mass ratio etc.). The detailed comparison between observa- 
tion and simulation (c.f. Figure W shows some systematic 
deviations which may be overcome by further improvements 
of the model. A possible cure for the above problems could 
be a numerical calculation of the field distortion in pres- 
ence of a highly conducting plasma and a more detailed 
investigation on the flow instabilities. However on account 
of the enormous Reynolds numbers of the flow a fully self- 
consistent numerical model seems not feasible to us. A fur- 
ther improvement for the model would be a more precise 
calculation of the radiation reprocessing rather that just as- 
suming limiting cases for the opacity, especially the velocity 
gradients in the stream could significantly alter the charac- 
teristics of the line radiation. However the overall agreement 
between observational data and the proposed model based 
on rather simplified physics seems remarkably good to us. 
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APPENDIX A: POWER SERIES FOR THE 
HORIZONTAL STREAM AT Li 

Vertical scale 

In their classical paper on gas dynamics in semi-detached 
binaries Lubow & Shu (1975) assume that thermal equilib- 
rium prevails in z-direction. In a later publication (1976) 
they give also a closed set of differential equations for the 
vertical structure of the stream: 

u 30 ^ + 2u X = (Al) 
as 

«s + -* + (S)„ = • < A2 > 

X(Q)=A, v(Q) = (A3) 
Therein is 1/-Jx the dispersion of the vertical Gaussian, 
v — ^f- the divergence of the z-velocity and ( ^f- ) the 
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second derivative of the Roche potential taken at the stream 
centre. 

Li is a singular point of these equations since energy 
conservation provides u s q = AiS + 0(s 2 ) and hence u s q is 
vanishing at initial point s = 0. The similarity to Euler's 
equation suggests the 'Ansatz' as a fractional power series 



X = A + 



k + c 



k+b 



(A4) 



for some b > 0, c > and sets {Ck} k and{iVfc}k with Co 7^ 
and iVo 7^ 0, which provides uniqueness of the solution. After 
substitution and use of Cauchy's theorem the s'-term in 
the x-equation and the s -term in the x-equation provide 
6 = c = 1, using the non- vanishing of Co and Nq. This 
proves that the solution we are looking for is even analytical 
at Li, which has not been obvious a priory. 

Knowing c and b we can compare the terms in s k and 
obtain to first order 
2A$ ZZS 



X = A + 



X( + 2A 
A? + 2A 



8 + 0(8") 
s + 0{s 2 ) 



(A5) 
(A6) 



where, with the white dwarf mass fraction fi, we have been 
using the definitions 



A := 



/' 



+ 



A? := 



i(A-2 + ^9A 2 -8a) 

/ <9 3 $ \ 
\d 2 zds) L , 



(A7) 
(A8) 
(A9) 



V9 2 n9s/ Ll 



(A16) 



APPENDIX B: MAGNETIC FOCUSING IN THE 
CURTAIN PLANE 

We now consider the magnetic focusing in the curtain plane. 
Be P a fixed point with co-latitude 9 on an arbitrary field 
line (r m , tp). On a neighbouring field line (r m + dr m , ip + dt/)) 
a point Q with the co-latitude 9 + d6 has a distance dl from 
P which is to first order in the locally Euclidean metric 



where 

r — 
dr — 



irf + (rddf + (rsmOd^f 



r m sin 9 

(r m + dr m ) sin 2 (9 + d9) — r — 
(2r m sin 9 cos 9)d9 + (sin 2 9)dr„ 



(Bl) 



(B2) 



(B3) 



The condition for a critical point is given by -^^(dl) 2 = 
and yields as a unique solution 



d9 = -2 



cos 9 sin 9 
r m (l + 3cos 2 9) 



dr„ 



(B4) 



which, if 9 > —d9, actually provides by substitution into 
Equation (Bl) the real minimum distance dl between P and 
any point Q(9 + d9) on a neighbouring field line: 



(dl) 2 = h 2 m (dr m ) 2 + h 2 b (dtp) 2 

h 2 _ 1 - cos 6 9 + 3 cos 4 9-3cos 2 9 

rm '~ 1 + 3 cos 2 9 

2 r m sin 2 6<(1 + cos 2 9) 

* := l + 3cos 2 



(B5) 
(B6) 

(B7) 



Horizontal scale 



Let 1 / ^7 be the horizontal dispersion of the Gaussian den- 
sity profile of the stream, /3 = the divergence of the 
n- velocity and K = 1/R the curvature of the ballistic tra- 
jectory. Then the differential equations for the horizontal 
stream width are, as derived by Lubow & Shu (1975): 

m s0 P- + 2/3 7 = 
ds 



(A10) 



Us ° + 132 ~ 7 = (2 ~ sKu ^ Kus ° - (13) - 4 ( An ) 



d 2 §- 



13(0) =0, 7(0) = A 2 — A + 3 =: A 2 



(A12) 



where we have already substituted the analytical integrals 
found by Lubow & Shu. The singularity of these equations 
is of the same kind as for the vertical structure, hence we 
are claiming again that 7 and (5 can be written as fractional 
power series. A somewhat similar but more extensive algebra 
shows then that the solution of the horizontal equations is 
also analytical. As first order solution we obtain 



7 = 



, 2 (12^oAr 
A 3 i — 



OA! 



(6A1 - 



2A 2 + A 2 
isjAi 



0(s 2 



2A| + A 2 



s + 0(s 2 



where Kg is the curvature at the Li-point and 



a; 



A? 



A + 2 



(A13) 



(A14) 



(A15) 



